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GLOBAL WEAK SOLUTIONS FOR THE THREE-DIMENSIONAL 
CHEMOTAXIS-NAVIER-STOKES SYSTEM WITH NONLINEAR DIFFUSION 

QINGSHAN ZHANG AND YUXIANG LI 


Abstract. We consider an initial-boundary value problem for the incompressible chemotaxis- 
Navier-Stokes equations generalizing the porous-medium-type diffusion model 

nt + u ■ Vn = An m — V • (nx(c)Vc), x € fi, t > 0, 

Ct + u ■ Vc = Ac — n/(c), x € 12, t > 0, 

Ut + k(u ■ V)w = Aw + VP + nV$, x € 12, f > 0, 

V ■ w = 0, ie!l, t > 0, 

in a bounded convex domain 12 C R 3 . It is proved that if m > |, ft G R, 0 < x U C 2 ([0,oo)), 
0 < / £ G * 1 ([0, oo)) with /(0) = 0 and <f> € W /1 ’°°(n), then for sufficiently smooth initial data 
(no,co, wo) the model possesses at least one global weak solution. 


1. Introduction 

Chemotaxis is the directed movement of living cells under the effects of chemical gradients. Aer¬ 
obic bacteria such as Bacillus subtilis often live in thin fluid layers near solid-air-water contact line, 
in which the swimming bacteria move towards higher concentration of oxygen according to mech¬ 
anism of chemotaxis and meanwhile the movement of fluid is under the influence of gravitational 
force generated by bacteria themselves. Both the oxygen concentration and bacteria density are 
transported by the fluid and diffuse through the fluid ( [5,14,20]). 

To model such biological processes, Tuval et al. [22] proposed the following model 

' nt + u ■ Vn = An — V • (n\(c)Vc), 
ct + u- Vc = Ac — n/(c), 
ut + k(u ■ V)n = An + VP + nV<f>, 

V ■ u = 0, 

in a domain fl C M^, where the scalar functions n = n(x, t ) and c = c(x, t ) denote bacterial density 
and the concentration of oxygen, respectively. The vector u = (ui(x,t),U 2 (x,t), ■ ■ ■ ,njv(x,t)) is 
the fluid velocity field and the associated pressure is represented by P = P(x,t). The function 
X is called the chemotactic sensitivity, / is the consumption rate oxygen by the bacteria and 
k € R measures the strength of nonlinear fluid convection. The given function 4> stands for the 
gravitational potential produced by the action of physical forces on the cell. 

The chemotaxis fluid system has been studied in the last few years and the main focus is on 
the solvability result. Under the assumption that x( c ) = X is a constant and / is monotonically 
increasing with /(0) = 0, Lorz [14] constructed local weak solutions in a bounded domain M. N 
(N = 2,3) with no-flux boundary condition and in M 2 in the case of inhomogeneous Dirichlet 
conditions for oxygen. In bounded convex domains C M 2 , Winkler [28] proved that the initial¬ 
boundary value problem for (1.1) possesses a unique global classical solution. In [30] the same 
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X € n, t > 0, 

illl, t > 0, 

x € n, t > o, 

t > 0 
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author showed that the global classical solutions obtained in [28] stabilize to the spatially uniform 
equilibrium (no,0,0) with no := py f n no(x)dx as t —>• oo. Zhang and Li [32] proved that such 
solution converges to the equilibrium (no, 0,0) exponentially in time. By deriving a new type of 
entropy-energy estimate, Jiang et al. [11] generalized the result of [30] to general bounded domains. 
For the well-posedness of the Cauchy problem to (1.1) in the whole space we refer the reader 
to [2,3,6,13,31,34], 

When the nonlinear convective term is ignored (n = 0 in (1.1)), which means the fluid motion is 
slow, and the model is simplified to the chemotaxis-Stokes equation. In this modified version, global 
weak solutions are constructed for the two-dimensional Cauchy problem [6]. In a bounded convex 
domain O C M 3 , the chemotaxis-Stokes system possesses at least one global weak solution [28]. 

The diffusion of bacteria may depend nonlinearly on their densities [9,19,20,23]. Introducing 
this into the model (1.1) leads to the chemotaxis-Navier-Stokes system with nonlinear diffusion [4] 


nt + u ■ Vn = V • (D(n)Vn) — V • (nx(c)Vc), 

iGSI, 

t > 0, 

ct + u ■ Vc = Ac — nf(c ), 

iGSI, 

t > 0, 

ut + k(u ■ V)u = A u + VP + nV4>, 

iGSI, 

t > 0, 

V • u = 0, 

x G Q, 

t > 0. 


Up to now, the main issue of investigation to (1.2) seems to concentrate on the chemotaxis-Stokes 
variant. Under the assumption D(ri) = mn m ~ 1 , Di Francesco et al. [4] proved that when m G 
(|,2] the chemotaxis-Stokes system admits a global-in-time solution for general initial data in the 
bounded domain SlcK". The same result holds in three-dimensional setting under the constraint 
rn G ( 7+ i 2 ^ , 2]. Tao and Winkler [19,20] extended the global existence result so as to cover the 
whole range m € (l,oo) in the bounded domain HcK 2 and m G (|,oo) in the bounded convex 
domain fl C M 3 . In [13], global existence of weak solution to the Cauchy problem of chemotaxis- 
Stokes system is established with m = ^ in fl = R 2 . Recently, Duan and Xiang [7] generalized the 
global existence result for all exponents m G [l,oo) . 

In contrast to the chemotaxis-Stokes system, very few results of global solvability are available 
for the full nonlinear chemotaxis-Navier-Stokes system. In the case U Q M 2 , global weak solutions 
are constructed by setting D(n ) = mn m_1 with m G [l,oo) [7]. For the three-dimensional initial¬ 
boundary value problem, the only result we are aware of is that when m > | the full system 
with nonlinear diffusion admits a global weak solution provided that <h G P 3 oc ((0, oo); Lj oc (Q)) with 
V<f> G -h/ 2 oc ((0, oo); L°°(Q)), and x and / are continuous differentiable satisfying x' > 0, / > 0 and 
/(0) =0 [24]. 

Recently, for sufficiently smooth initial data (no, co, uo), Winkler [29] established global weak 
solutions of (1.1) in bounded convex domains H C R 3 under the assumptions x £ C^QO, oo)), 
/ G C' 1 ([0, oo)) with /(0) = 0 and <h G VF 1,00 (fI). Motivated by the work of [29], our purpose of the 
present paper is to consider the full chemotaxis-Navier-Stokes system with nonlinear diffusion. In 
order to formulate our result, we specify the precise mathematical setting: we shall subsequently 
consider (1.2) along with boundary conditions 

f)Tl fjr 

D(n)— = — = 0, u = 0, x G d£l, t > 0 (1.3) 

on on 


and the initial conditions 


n(x, 0) = no(x), c(x, 0) = cq(x), u(x, 0) = uq(x), iGSI 


(1.4) 
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in a bounded convex domains Q C M ' 3 with smooth boundary, where we assume 
no £ LlogL(fi) is positive, 

Co € L°°(fi) is nonnegative and such that y/co € IF 1 , 2 (n), ( 1 . 5 ) 

Uo £ L 2 (fl). 

With respect to the parameter function in (1.2), we shall suppose throughout the paper that 


B(*)eC(( 0 , 00 )) for some 7 > 0 , 
D lS m ~ l < D(s) < Dis™- 1 for all s > 0 


with m > | and D 2 > D\ > 0, and that 


X € C 2 ([ 0 , 00 )), X > 0 in [ 0 , 00 ), 

/ € C' 1 ([ 0 , 00 )), /( 0 ) = 0 , / > 0 in ( 0 , 00 ), 

4> <E W ,1 ’ 0O (fi). 


Moreover, we shall require the further technical assumptions 

— ^ > 0 , on [ 0 , 00 ) 


< 0 , on [ 0 , 00 ) 


and 


(x-f)' > 0 , on [ 0 , 00 ). 


( 1 . 6 ) 

(1.7) 


( 1 . 8 ) 


(1.9) 

( 1 . 10 ) 

( 1 . 11 ) 


Our main result reads as follows. 


Theorem 1.1. Let Ll C M 3 be a bounded convex domain with smooth boundary and k £ M. Suppose 
that the assumptions (1.5)-(1.11) hold. Then there exists at least one global weak solution (in the 
sense of Definition 6.1 below) of (1.2)-(l-4) such that 

£ TL([°>o°); W l ’ 2 (Sl)) and c^ € Tf oc ([0, 00 ); W M (^))- 

Remark 1.1. (i) If the diffusion function D{u) = 1 in (1.2), this is consistent with the result 
of [29]. 

(ii) Theorem 1.1 shows that the model (1.2)-(1.4) possesses global weak solution even when the 
diffusion effect is rather mild. However, we have to leave open here whether the lower bound of 
diffusion exponent m = | is optimal to guarantee global weak solvability. 

The rest of this paper is organized as follows. In Section 2, we introduce a family of regularized 
problems and give some preliminary properties. Based on an energy-type inequality, a priori esti¬ 
mates are given in Section 3. Section 4 is devoted to showing the global existence of the regularized 
problems. In Section 5, we further establish some e-independent estimates. Finally, we give the 
proof of the main result in Section 6 . 

Notations. Throughout the paper, for any vectors v £ M' 3 and w £ M' 3 , we denote by v <g) w the 
matrix ^. 3 x 3 with ay = ViWj for i.j £ {1,2,3}. We set LlogL(H) is the standard Orlicz space and 
L 2 (H) := {<£> £ L 2 (0)|V • tp = 0 } denotes the Hilbert space of all solenoidal vector in L 2 (Q). As 
usual V denotes the Helmholtz projection in L 2 (H). We write Wq'^(D) := Wo’ 2 (fi) H T 2 (H) and 
:= Cfi°(Q) n L 2 (H). We represent A as the realization of Stokes operator —VA in L 2 (H) 
with domain D(A) := W 2 ’ 2 {fii) n Wq ^(VL). Also n(-,t), c(-,t) and u(-,t) will be denoted sometimes 
by n(t), c(t) and u(t). 




4 


ZHANG AND LI 


2. Regularized problem 

Our intention is to construct a global weak solution as the limit of smooth solutions of appro¬ 
priately regularized problems. According to the idea from [29] (see also [20]), let us first consider 
the approximate problems 

n £t + u £ • Vn £ = V • (D £ (n £ )\7n £ ) - V • {n £ F' e (n £ )x(c £ )S7c £ ) 
c £ t + u £ ■ Vc £ = Ac £ - F £ (n £ )f(c £ ), 
u £t + k(Y £ u £ ■ V)u £ = A u £ + VP £ + n £ V<f>, 

V ■ u £ = 0, 

§H£. — ®£e. — f) v — 0 
du ~ du ~ U > u e ~ U ’ 

n £ (x, 0) = n 0£ (x), c £ (x, 0) = c 0£ (x), u £ (x,0) = u 0£ (x), 

< 

for e G (0,1), where the approximate initial data no £ > 0, co £ > 0 and uq £ satisfy 

' n 0e e 

< fn n 0e = fn n 0 , ( 2 . 2 ) 

no £ —>• no, £ —>• 0 in L log L(R), 

' v^eW), 

< ll c 0 e||z,°°(Q) < I|co||l°°(D), (2.3) 

y/coe -4 y/coi £ —>• 0 a.e. in R and in IY 1 , 2 (fl), 

and 

' UOe G Cg^(n), 

< II^OeIIlRU) = ll«o||L 2 (n), (2.4) 

woe —>• U 0 ) e —>• 0 in T 2 (fl). 

The approximate functions in (2.1) can be chosen as 

D £ (s ) := D(s + e), for all s > 0, (2-5) 

F £ (s) := - ln(l + es), for all s > 0, 
and the standard Yosida approximate Y e ( [17]) is defined by 

Y £ v := (1 + eA)~ l v , for all v G L 2 (R). 

It is easy to verify our choice of F e above guarantees that for each £ £ (0,1) 

0 < F'(s) =- < 1, for all s > 0, (2.6) 

_ £W 1 +£S V ’ 

sF'(s ) = —-— < -, for all s > 0, (2.7) 

1 + es £ 

0 < F £ (s) < s, for all s > 0, (2.8) 

and 

F e (s) —> 1, T £ (s) —> 1, as e —>• 0 for all s > 0 . 

The first lemma concerns the local solvability of the approximate problems (2.1). The proof is 
based on well-established methods involving the Schauder fixed point theorem, the standard regu¬ 
larity theory of parabolic equation and the Stokes system (for details see [20,28,29], for instance). 
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Lemma 2.1. For any e G (0,1), there exist a maximal existence time T maXi£ € (0, oo] and deter¬ 
mined functions n £ > 0, c £ > 0 and u £ 

n £ € C° (fi x [0,T maXi£ )) n C 2,1 (Cl x (0,T maX]£ )) , 

c £ € C° (Cl x [0, r maX;£ )) n C 2 ’ 1 (Cl x (0, T maX;£ )) and 

u £ eC°(Clx [0,T max>e )) n C 2 ’ 1 (fi x (0, T maXj£ )) 

such that (n £ ,c £ ,u £ ) is a classical solution of (2.1) in Ll x (0, T maXj£ ). Moreover, if T maX]£ < oo, 
then 

ll^e(') IIl°°(Q) “I" ll c £ (') Iliy 1 '9(f2) “I" ||^e(‘j ||d(A“) t ^ t ^max,E 
/or all q > 3 and a > 

The following estimates of n £ and c £ are basic but important in the proof of our result. 


Lemma 2.2. For each e G (0,1), we have 

/ n £ (-,t) = n 0 for all t G (0,T maXj£ ) 
Jn Jn 


and 


\ce (’)i)||L°°(f 2 ) — I|co||l oo (o) — • M in n x (0,T maXi£ ). 


(2.9) 

( 2 . 10 ) 


Proof. Integrating the first equation in (2.1) and using (2.2), we obtain (2.9). Since / > 0 by 
our assumption (1.8) and F e > 0 by (2.8), an application of the maximum principle to the second 
equation in (2.1) gives (2.10). □ 


3. An energy-type inequality 


In this section, we shall utilize an energy inequality associated with the first two equations in (2.1) 
to establish a priori estimates. The inequality is frequently used in the literature (see [6,20,29,30], 
for example) and it also will play an important role in our proof. 

Lemma 3.1. Let (1.6)-(1.11) hold. There exists K > 1 such that for any s G (0,1), the solution 
of (2.1) satisfies 


d 


dt 


1 


-77 i / n £ lnn £ + - / |VT(c £ )| > + — 


in 


in 


1 


D £ (n £ ) 


K [Jn n E 


|Vn e | 2 + / ■ 
Jn 


| D 2 c 


+ 


|Vc £ 


n 4 


< K I | V| 2 for all t G (0,T maXi£ ), 
Jn 

where T(s) := ff If— with g(s ) := ; 44. 
v ’ Jl -/gM v ' x(s) 


(3.1) 


Proof. The proof is based on the first two equations in (2.1) and integration by parts and detailed 
computations can be found in [29, Lemmas 3.1-3.4]. □ 

Based on Lemma 3.1, we can modify the above energy-type inequality (3.1) to contain all com¬ 
ponents of n £ , c £ and u £ . 

Lemma 3.2. Let T be as given by Lemma 3.1 and suppose that (1.6)-(1.11) hold. Then for any 
(0,1), there exists C > 0 such that 


s{/A‘ ln "' + U |v * fe)|2 Wj“ !|2 } 


+ 2K' Dl I + l V ^l + 


\D 2 c £ 


+ 


in 


[ I Vc £ 

<n 4 


I Vu £ 


<C for all t G (0,T maX!£ ), 

where D\ and K are constants provided by (1.7) and Lemma 3.1, respectively. 


(3.2) 
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Proof. Multiplying both sides of the third equation in (2.1) by u £ and integrating by parts over U, 
we have 

f Kl 2 + / |Vu £ | 2 = -k f (Y £ u £ -\7)u £ -u £ + f n £ V4> • u £ for all f € (0,T max , £ ). 

2 at Jq Jn Jq Jn 

Since V • u £ = 0 implies V • Y £ u £ = 0, we thereby obtain 

“ I \u £ \ 2 + [ |Vu £ | 2 = f n £ V<h ■ u £ for all t £ (0,T maXi£ ). 

2 at jfi Jq Jq 

Substituting this into (3.1), we get 

f t {f a n^n c + \f a \^ { c.)? + K jwf} 


+ K 


Iff D £ (n £ ) 


|Vn £ | 2 + 


\D 2 c t 


2 „ |2 
£\ 


+ 


In n e Jn c e 

< 2K f n £ Vd> • u £ for all t £ (0,T maX)£ ). 
Jn 

Using (1.7) and (2.5), we have for each e € (0,1) 
f D £ {n £ ) 


f |Vc £ |^ 

In 4 


+ K / |V« e | : 


In n e 


IVn e | 2 > D l [ (n £ + e) m_2 |V?r e | 2 for all t £ (0,T max , £ ). 
Jn 


(3.3) 


(3.4) 


By (1.8), Hoder’s inequality and the embedding PU 1,2 (U) L 6 (U) for n = 3, we can find C\ > 0 

such that for each e £ (0,1) 

2K Sn^^ Ue ~ + e|l L | (n )||^||L 6 (Q) 

< Ci\\n £ + e|| i 6 (n) ||Vu £ || L 2 (r2 ) for all t £ (0,T maX)£ ). 


Note that 

||(n £ + e)T|| i i (n) < (||reo|| L i(r2) + |U|) 2 for all t £ (0,T maX)£ ) 
by (2.2) and (2.9). It follows from the Gagliardo-Nirenberg inequality [25] that 


(3.5) 


l n ' + e llzJ (n ) 


, x EL 

(n £ + £) 2 


12 

L^{Q) 


< C 2 


m 

V (n e + e) 2 


L 2 (Q) 


(n £ + e) 2 


5m — 2 
m(3m-1) 

2 i 

L rn ($~2) 


/ \ Ill 

(n £ + e) 2 


L m (Q) 


< C 3 [ ( ^(n £ + e) m - 2 |Vn £ | : 


2(3m—1) 


+ 1 for all t £ (0, T maXj£ ) 


with C*2 > 0 and U3 > 0. Since 3m 1 _ 1 < 1 by our assumption m > | , Young’s inequality yields 
C4 > 0 such that 


2 K / n £ V4>-u £ < CrC 3 / (n £ + e) m ' 2 |Vn £ | 


/n 


2(3m-l) 


1 

K 

Di 


Y —Cl\C'i ( ( I (n £ + e) m 2 |Vn £ |^J + 1J + 

< ^ / (n £ + er- 2 |Vn £ | 2 + ^- f |V« £ | 2 + C- 4 


+ lj • ||V« e || L 2 (n) 
K 


|Vu £ |" 


for all t £ (0,T maX)£ ). Inequality (3.2) then follows by combining (3.3), (3.4) and (3.6). 


(3.6) 

□ 
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We can now use Lemma 3.2 to establish a priori estimates of the solution of (2.1). 

Lemma 3.3. Let 'I' and K be as given by Lemma 3.1, and assume that the requirements of Lemma 
3.2 are satisfied. Then there exists C > 0 such that for any e € (0,1) we have 


n e \nn e + \ f |VT(c £ )| 2 + K [ \u £ \ 2 < C for all t £ (0,T maX!£ ) 
> J Jo Jo 


(3.7) 


and 


Dl / 

Jo Jo. 

for all T € (0, T maX;£ ) 
Proof. Set 


(n £ + e) m_2 |Vn £ | 2 + f f P [ |Vu £ | 2 < C(T + 1) (3.8) 

Jo Jn c e Jo Jn c e Jo Jn 


Veifi) \= [ n £ lnn £ + ^ [ |W(c e )| 2 + if [ \u £ \ 2 for all t G (0, T max>e ) 
Jo. z Jo Jo 


(3.9) 


and 


D 2 c' 2 


h £ (t ) := D i / (n £ + e) m |Vn £ | 2 + 

Jo Jo Cc 


+ 


f |Vc £ | 4 

In c a 


+ / |Vu e | 2 for ah t € (0, T maXi£ ). 


Then (3.2) implies 


y £ (t) + — h £ (t) < C for ah t € (0, T maX)£ ). 


(3.10) 


In order to introduce dissipative term in (3.10), we show that y £ {t) is dominated by h e (t). Now 
using the inequality 

3 


z In 2 i < 


3m — 1 


-z m+ 3 for all z > 0 


with m > we can find positive constants C±, C 2 and C 3 fulhlling for each e € (0,1) 


/ 


n F In n F < - 

n 3m - 1 


(n £ + e) m+ l 


< 


3m — 1 

3Ci 
3m — 1 


m £ + e 2 


2(3m+2) 

3m 

2(3ra+2) 

L——(f2) 


\ m 

V(n £ + e) 2 


3m 

3m+2 


< C 2 


. m 

V (n £ + e) 2 


L 2 (n) 


L 2 (Q) 

/ \ m 

(n £ + e) 2 


/ \ 111 

(n e + £) 2 


2 

3m+2 

2 

L rn 


4 

3m 

2 1 

L rn 


(n £ + e) 2 
(n £ + e) 2 


L m (Q) 


2(3m+2) 

3m 


2(3m+2) 
3 m 


L m (Q) 


< C 3 ( ^(n £ + e) m - 2 |Vn £ | 2 + l ) for allt £ (0, T max , £ ) 


(3.11) 


by the Gagliardo-Nirenberg inequality and (3.5). According to (1.8), we have 
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Hence the mean value theorem yields g(s) > min tg[o,m]9'{ t ) s =: Ms. We now apply Young’s 
inequality and ( 2 . 10 ) to obtain 


U ™ 2 = \L 


< 


|Vc £ 

g(c £ 


1 f |Vc £ 


< - 


In c e 

1 f I Vc £ 


+ 


+ 


4 Jn 9 2 (ce) 


M\n\ 

4 M 2 


4 Jn ct 

From the Poincare inequality, we have C 4 > 0 such that 

|2 


for all i € (0, T n 


/L [ \u £ \ 2 < C 4 f |Vu £ | 2 for all t G (0,T maXj£ ). 

JO IS) 

It follows easily from (3.11)-(3.13) that 

Ve(t) < C 5 h e (t) + C 6 for all t G (0,T maX;£ ) 

with C 5 := max C 4 j and C§ := C 3 + This, along with (3.10), yields 

Veit) + 4 ^r2/e(*) + (*) < CV := C + for all t € (0,T maXj£ ). 

Noting that h £ (t) > 0, a standard ODE comparison argument implies 


(3.12) 


(3.13) 


Veit) < max l sup y e (0), 4KC 5 C 7 \ for all t G (0,T maXj£ ). 
Ug(o,i) I 


(3.14) 


In view of (2.2)-(2.4) and [29, Lemma 3.7], we obtain (3.7). On the other hand, since zlnz > -- for 

.M 

e 


all z > 0, we have y £ (t) > — —. Therefore, a time integration of (3.14) directly leads to (3.8). □ 


4. Global existence for the regularized problem (2.1) 

With Lemma 3.3 at hand, we are now in the position to show the solution of approximate problem 
( 2 . 1 ) is actually global in time. 

Lemma 4.1. For each e G (0,1), the solutions of (2.1) are global in time. 

Proof. In this section, we shall denote by C various positive constants which may vary from step to 
step and which possibly depend on e. Assume for contradiction that T max>£ < 00 for some e G (0,1). 
By Lemma 3.3, we know that 


In 


\u £ \ 2 < C for all t G (0, T maXi£ ) 


(4.1) 


and 


rF nax,e r r'l 

/ / |Vc £ | 4 = 

Jo Jn Jo 


Me? 

In 4 


f*T. max,e P 

< M 3 / 

Jo Jn 


|Vc £ 


„ r 3 

n G 


< C for all t G (0,T maX)£ ). 


(4.2) 
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Multiplying the first equation in (2.1) by p(n £ +e) p 1 withp G [m+1, 2(m+l)] and using integration 
by parts we obtain 

4 [ (n £ + e) p +p{p- 1) f (n £ + e) p ^ 2 D £ (n £ )\Vn £ \ 2 


= P(j> - !) / fae + e) p 2 n £ F' e (n £ )x{c £ )Vc £ ■ Vn £ 

Jn 

for all t G (0,T maX)£ ). We deduce from (1.7), (2.7) and Young’s inequality that 
4 f (n £ + e) p + p(p - 1)D 1 f (n £ + e) m+p “ 3 | Vn £ | 2 
<p(p —!)_ max f (n e + e) p_ 2 Vc £ ■ Vn £ 

£ s£[0,M] 


(4.3) 




< 


p(p-l)I>i [ (n £ + e) m+p - 3 \Vn £ \ 2 + [ (n £ + £ ) 2 (p—D + C [ |Vc £ 


for all t G (0,T maX)£ ). Since 2(p — m — 1) < p for p < 2(m + 1), applying Young’s inequality again, 
we obtain 

4 /4h + e) p < [ (n £ + e) p + C I | Vc £ | 4 + C for all t G (0, T max , £ ). 

Integrating this yields 


f (n £ + £) p < C for all t G (0, T maXi£ ), 
Jn 


(4.4) 


where p G [l,2(m + 1)]. We now use the idea from [29] to obtain the boundedness of u £ . From 
(4.1), we get 

||Y £ « £ (t)|| L oo (a) = ||(1 +£A)~ 1 u £ (t)\\ L o 0{n) 

< C\\u £ (t)\\ L 2 m 

< c for all t G (0, T maXj£ ) (4.5) 

due to the embedding 11(1 + eA) L°°(Q). We apply the Helmholtz projection V to the third 
equation in ( 2 . 1 ), test the resulting identity by Au £ and integrate by parts over H to have 

|^w | 2 + / l^e ] 2 = [ VH £ ■ Au e for all t G (0, T maXi£ ) 

Jn Jn 


2 dt Jn 

with H e (x,t ) := n £ V<h — k(Y £ u £ ■ V)u £ , where we have used 


[ (j) ■ A(j> = [ \A l 2 (t)\ 2 = [ |V0 | 2 for all G D(A). 
Jn Jn Jn 


(4.6) 


Applying Young’s inequality, \\V(/)\\ L 2 ^ < ||^||l 2 ( q) for all (j) G L 2 (Q) ( [17, Lemma 2.5.2]), (4.5) 
and ( 1 . 8 ), we can estimate 

[ VH £ ■ Au £ < [ \Au e \ 2 + \ [ \VH e \ 2 

Jn Jn 4 Jn 

< \ \Au e \ 2 + ^- [ \{Y £ u £ -Y)u £ \ 2 + \ [ |n £ V$| 2 

Jn 1 Jn 1 Jn 

< [ \Au £ \ 2 + c( f |Vu £ | 2 + [ n 2 ] for all t G (0,T maXi£ ). 

Jn \Jn Jn J 


Hence we get 


1 d_ 

2 dt 


[ \A2u £ \ 2 <c( [ |Vu £ | 2 + [ n 2 ] for alU G (0,T max>£ ). 

Jn \Jn Jn J 
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This, along with (4.4) and (4.6), gives 

[ |Vu £ | 2 < C for all t G (0,T maXi£ ). 

Jn 

We thereby obtain 

\\VH £ (t)\\ L 2 {n) <C for all t G (0, T maX)£ ). (4.7) 

Applying the fractional power A a with a G (|, 1) to both sides of the variation-of-constants formula 

u e {t) = e~ tA u 0£ + f e^ t ~ s ' )A VH £ (s)ds for all t G (0,T maXj£ ), 

Jo 

using the well-known smoothing estimate of the Stokes semigroup ( [ 8 ]) and (4.7), we have 

||A a 'u e (t)|| £ 2 (Q) < || A a e tA uo £ \\L^(n) + f \\A a e ^ b ' lA 'PH £ (s)\\ L 2 ^ds 

Jo 

< Ct a \\u 0 e \\ L 2 (ty + C f (t — s) a \\VH e {s)\\ L 2 (^ds 

Jo 

< C for all t G (r, T maX;£ ) 

with any r G (0,T max £ ). In view of the embedding D(A a ) L°°(Q) asserted by our chioce of a 
( [17, Lemma 2.4.3]), we deduce 

IK(*)IU°°(fi) < C\\A a u e (t) || L 2 ( n) < C for all t G (r,T maXj£ ). (4.8) 

Let r := min{2(m + 1), 4}, then r > 3 due to m > |. Employing V to both sides of the variation- 
of-constants formula for c £ 

c £ (f) = e^ - i )A c £ (^) - J e (t_s)A (F e (n e )f(ce) + u £ • Vc £ ) (s)ds for all t G (^, T maX)£ ), 

recalling the standard smoothing estimates of Neumann heat semigroup ( [27, Lemma 1.3], see 
also [15]), (2.8), (2.10), (4.8), (4.4), (4.2) and the Holder inequality we obtain 

rt 


m 1 c £l)IIz/'(G) ^ 


Ve^ )A c e (-) 


L r (Cl) 


Ve o- S )A (F£ ( ne ) / ( Ce ) + Ue . Vce ) ( s ; 


L r (fl) 


ds 




c £ (~) . + C (*“ s ) ’ ||(iW(ri £ )/(c e ))(s )||^(^1 

2 ZA(f2) Jr 


+ c j r (t-s) 3 H(u £ - Vc £ )(s)|| Lr(n) ds 


< Ct 2 +C I (t-s) 2\\n £ (s)\\ L r( Q) ds + C I (t - s) 2 ||Vc £ (s)|| iI . (Q) ds 


(«)' 


< 


Cr 2 +cj(t — s ) zds + C 1 J (t — s ) 2 ||Vc £ (s) 11^4^) 


ds 


< Cr ”2 


^ 1 / r fmax,£ 2 \ 4 / /*^max,e /* 

+ ( t- s )—3dsj U i v c £ (s)| 


< Ct 2 + CTmax^ + CTmax., £ 


< C for all f G (r,T maX)£ ). (4.9) 

! next rewrite the variation-of-constants formula for c £ in the form 

c £ (f) = e t(A_1) c 0 , e + [ e ^ _s)(A_1) (c £ - F £ (n £ )f(c £ ) - u £ • Vc £ ) (s)ds for all t G (0, T maXi£ ), 

Jo 
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Picking 0 G (^ + ^,1), then the domain of the fractional power D((—A+l) y ) VP 1 , 0 O (0) ( [26,33]). 

Hence, we obtain by virtue of L p -L q estimates associated heat semigroup ( [26]), (2.10), (2.8), (1.8), 
(4.4), (4.8) and (4.9) 


l c e( i )llvE 1 '°°(n) < C (—A + l) e c £ (t) 


L r (tt) 


< cr e e- vt ||co, e | 


L r (f2) 


+c 


/ (t-s) 
Jo 


_ 0 e — s ) 


(c £ - F £ {n £ )f{c £ ) - u £ ■ Vc £ ) (s)|| Lr(n) ds 


< CT~ e + C j\t - s)- e e- v{ t- s) ds + C J\t - s)- e e- u{ t~ s ) ||n e («)|| Lr(n) ds 




m+p— 3 1 
L e 


n; 1 "'' “|Vn £ | 


< CT~ e + CT{l-Q) 

< C for all t € (t, T maX;£ ) (4.10) 

with some p > 0, where T(-) is the Gamma function. We may then apply (4.3) once more and 
Young’s inequality to deduce that 

if a nt + P(r- l)A/ o 

< C [ n p ~ 2 Vc £ ■ Vn £ 

Jn 

< C f nr 2 |Vn e | 

Jn 

<p(p-l)D 1 [ n ™+P- 3 |Vn e | 2 + [ n p - m - l + C 

<p(p-l)Di [ n™ +p ~ 3 \Vn £ \ 2 + [ n p + C for all f € (r, T maXi£ ). 

Jo, Jci 

Therefore, integrating with respect to t, we obtain 

I n p <C for all t G (r, T maX]£ ) 

Jn 

with any p > 1. Upon an application of the well-known Moser-Alikakos iteration procedure ( [1,18]), 
we see that 

IK(*)IU°°(n) < C for all t €. (t 5 ^max,e ). (4.11) 

In view of (4.8) (4.11), we apply Lemma 2.1 to reach a contradiction. □ 

5. Further e-independent estimates for (2.1) 

In order to pass to limits in (2.1) with safety, we need some more e-independent estimates for 
the solution. 

Lemma 5.1. Suppose that (1.6)-(1.11) hold. There exists C > 0 such that for all e € (0,1), the 
solutions of (2.1) satisfy 

r T 


II 


Vn £ : 


< C(T + 1), for all T > 0, 


(5.1) 


/ 7 < 


n £ + e) p <C(T + 1), l<p< 


3 m + 2 


3 


for all T > 0, 


(5.2) 
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ff a 


D e (n e ) 


I o Jn n e 

rT 


|Vn £ | 2 < C{T+l) for all T> 0, 


// 

Jo Jn 


3m-)-2 


(D £ (n £ )Vn £ ) 3m + 1 < C(T + 1) for all T > 0 


and 


ff 

Jo Jn 


u £ | 3 < C(T + 1) for all T > 0. 


Moreover, if g < m < 2, t/ien we have 

f 3m+2 

/ / |Vn £ |^“ < C(T + 1), for all T> 0. 

Jo jn 

Proof. From Lemma 3.3 we know that there exists Ci > 0 such that 

rT 


[ [ (n e + e) m_ 2 |Vn £ | 2 <C\(T+ 1) for all T > 0. 

Jo Jn 


(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


Then, (5.1) is a direct consequence of (5.7). Due to the fact that D is bounded, we only need to 
prove (5.2) with p = 3m + 2 . We employ the Gagliardo-Nirenberg inequality to find C 2 > 0 and 
C 3 > 0 fulfilling 


[ T [ {n £ + e) 
Jo Jn 


3m+2 

3 = 


/ x m 

(n £ + e) 2 


2(3m+2) 

3m 

2(3m+2) 

L 3 In ( Q ) 


< c 2 ( 

\ m 

V(n e + e) 2 

3m 

3m+2 

L 2 ( n ) 

/ x 2Z1 

(n £ + e) 2 

2 

3m+2 

2 1 

l "i (r2) 

/ X 222. 

(n £ + e) 2 

II 

S 

. m 

V (n £ + e) 2 

2 

L 2 (fi) 

m 

(n e + e) 2 

4 

3m 

2 1 

L ra (f2) 

m 

(n £ + e) 2 


2(3m+2) 

3m 


L m (f2) 

2(3m+2) ' 
3 m 

LTn (Q) 


< C 3 (r + 1) for all T > 0. 
Recalling the proof in Lemma 3.2, we obtain 


m+ w<L 


D £ (n £ ) 


|Vn e | 2 < K for all T > 0, 


In n e 

where y e and K are provided by (3.9) and Lemma 3.1, respectively. Integrating this in time over 
(0, T) yields (5.3). By Holder’s inequality, (5.2) and (5.3), we can find C 4 > 0 such that 

3m+2 3m 

[ [ (D £ (n £ )Vn £ )“ < ( [ [ Vn £ | 2 ) ^ ( [ [ ( D £ (n £ )n £ ) 

Jo Jn \J 0 Jn n e J \J 0 Jn 

D £ (n £ ) 


3m+2 
3 m 


/a 

3m-|-2 / pT 


6m+2 


< D % m+2 ( f f CK '^' \ Vn £ 

J0 JG 


3m + 2 

6m+2 / /*-* 


0 


/O JG n £ 

< C 4 (T + 1) for all T > 0. 

Using the Gagliardo-Nirenberg inequality and Lemma 3.3 again we have 

rT 


f 3m+2 

/ (n £ + e) 3 
Jn 


6m+2 


fi 


10 

3 

. 10 , x 
L“3"(f2) 


[ M*)|| 

J0 

< C 5 J^ (||Vu eW n L 2 (Q) 

< C 6 (T+1) for all T > 0 


£\ u )\\ 2 TMn\ ‘ ll^e(^) II 3 


L 2 (n) 


— \ 

+ ll u e( i )lll 3 2(Q) ) 
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with C 5 > 0 and Cq > 0. Finally, we prove (5.6). Since | < m < 2, applying (5.7) and Young’s 
inequality we get C 7 > 0 and C 8 > 0 such that 


[ T [ I Vn £ 
Jo Jo 


3m-|-2 


p (m —2)(3m+2) 3m+2 (2 — m)(3m+2) 

= n £ 8 |Vn £ | 4 n £ 

< C 7 ([ T [ n?~ 2 \Vn £ \ 2 + f [ ) 

\J 0 Jfi 7o Jn J 

< C 8 (T + 1 ), for all T > 0 . 


This completes the proof. 


□ 


We derive an L p -bound for n e + e and a estimate of space-time integral / Q T | V(n £ + e ) r | 2 as a 
supplement to the regularity property concerning n e in the case m > 2 . 


Lemma 5.2. Let m > For all e € (0,1), there exists C(T ) > 0 and C > 0 suc/i that 


f (n £ + e) p < C'(T) /or all t > 0. 
Jn 


with 1 < p < 9 (m — 1) and 


with y < r < 5 (m — 1). 


|V(n £ + ef | 2 < C(T + 1) for all T> 0. 


(5.8) 


(5.9) 


Proof. It is based on a bootstrap argument ( [20], see also [10,21]). We multiply the first equation 
in ( 2 . 1 ) by p(n £ + e ) p_1 to deduce that 


4 [ (n £ + e) p + p(p - 1 ) f (n £ + e) p 2 D £ (n £ )|Vu 

= p(p - !) [ n £ (n £ + e) p_ 2 F £ (n £ )x(c £ )Vc £ • Vre ; 
Jn 


(5.10) 


for all f > 0. However, unlike the proof of Lemma 4.1 we deal with n £ F' e {n £ ) together, because our 
goal is to get an e-independent bound (5.8). More precisely, from (5.10) and (2.6) we have 


d_ 

dt 


[ (n £ + s) p + p(p - l)Di f (n £ + e) m+p 3 | Vn £ 
Jn Jn 

8 p(p- 1 ) max x( s ) / {n £ + eY^ 1 ' 
se[o,M]' Jq 


Vc £ • Vn £ 


(5.11) 


for all t > 0. Applying the Holder and Young inequalities in the right-hand side of (5.11), we have 
C 1 > 0 such that 


d_ 

dt 


'n 


(n e + ef + p(p - 1)D! / (n £ + e) m+p_3 |Vr 

Jn 


f \ f \ / \ rn + P— 3 

^ P\P — 1) max x{ s ) \ n e + £ ) 2 Vn £ 
sG[0,M] ' 


L 2 (n) 


p-m+1 

(n £ + e) 2 


L 4 (Q) 


|Vc £ || L 4 (n) 


< 


P(P ~ 1)-Di 


f (n £ + e) m+p- 3 |Vn e | 2 + Ci (n £ + e) 
Jn 


p— m+1 
2 


L 4 (n) 


l^ 7c elli 4 (n) (5.12) 


for all f > 0 . Assume that 


/ (n £ + e) Pi < C(T) 
Jn 


for all f > 0 























14 


ZHANG AND LI 


holds with some pi > 1 (this is true for p\ := 1 by (2.9) and e < 1). The Gagliardo-Nirenberg 
inequality gives C 2 > 0 such that 


, p-m+1 

2 

. p+m-l 

(n £ + e) 2 

L 4 (f2) 

(n £ + e) 2 


2(p—m+1) 
p+m — 1 
4(p — m+1) 

L P+" 1 - 1 (fi) 


< Cal 


+ 


_ , . P+m-l 

V(n £ + e) 2 


a 

L 2 (ft) 


p+m—1 


(n £ + e) 2 


1—a 

£ p+m-l (Q) 


p+m-l 

(n £ + e) 2 


2(p-m + l) 
p+m —1 


2 gj 

ip+m _! (Q) 


for all t > 0 , where 


If 


1 — Pi 

^TTT e( 0 ,l). 

1 3(p+rn-l) 


2(p — 771 + 1) 
p + m — 1 


a = 1 , 


which is equivalent to p = 3(m — 1) + |pj, and therefore by Young’s inequality 


C 1 


p—m+1 


(n £ + £■) 2 


L 4 (f2) 


|Vc, 


IIl 4 (g) 


<C 3 

p(p- l)T»i 


_ . . P+m-l 

V (n e + e) 2 


< 


L+G) + 1 ) l|VCell ^( n ) 

(„ e + e)m+P -3| Vn£ |2 + C 4 ||Vc £ ||t 4(a) + C 5 

4 JQ 


for all t > 0 with certain positive constants C 3 , C 4 and C 5 . Substituting this into (5.12), we obtain 

I J (,n, + 0” + h£_i+l („, + e) ™+»-3| Vn € \ 2 < C 4 ||Vc«|||, (n) + C 5 (5.13) 

for all t > 0. By integration, we finally get 

[ ( n £ + e) p < C(T) for all t > 0. 

Jn 

with p = 3 (m — 1) + |pj. By this iterative procedure, there exists a sequence {p*} such that 

, , 2 
Pi+l = 3(771 - 1) + ~Pi- 

It is easy to check that the sequence {p*} is increasing and pi —>• 9(m — 1) as i —>• oo. Therefore, 
we can reach any p < 9(m — 1) by finite steps and (5.8) is thereby proved. Another integration of 
(5.13) yields 


[ [ V(n £ + e)^ 

Jo Jn 


< 


Cq ( [ (jio £ + 1 Y + [ 
\Jn Jo 


IV c e 1 4 + 1 


in 


< C 7 {T + 1) for all T > 0 
with C 6 >0 and C 7 > 0 . This proves (5.9). 


□ 


In order to derive strong compactness properties, we also need some estimates concerning the 
time derivative of the solution. 
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Lemma 5.3. Let 7 := max{l, }. There exists C > 0 such that for all £ € (0,1) we have 


with some q > 3, and 


r 

d 

—ril 

Jo 

dt £ 

r T 

d _ 


dt ^ 

r 

d 

L 

ot ^ 


dt < C(T + 1) 

/or all T > 0 

(5.14) 

(tE 2 '9(n))* 



5 

dt < C(T + 1) 

/or all T > 0 

(5.15) 




5 

5 df<C(T+l) 

... 

/or all T > 0. 

(5.16) 


as well as 


Proof. Multiplying the first equation in (2.1) by yrzj -1 ^ with (p € C°°(fl) and integrating by parts, 
we obtain 


(n])tV 


- 7 ( 7 - 1 ) / D £ (n £ )nf 2 |Vn £ |V-7 / D £ (n £ )nf l Vn £ -Vy 


+7(7 - 1 ) [ nj 1 F £ (n £ )x{c £ )Vc E • Vn £ (p 
Jn 

+7 [ nJF £ (n £ )x(c £ )Vc £ ■ Vip+ [ nfu £ ■ Vp 
Jn Jn 

< ("7(7-1) [ 71 £ (n £ )n/ _ 2 |Vn £ | 2 + 7/ Tl^n^rc^lVr+l 

V Jn Jn 

+7(7-1) f F £ (n £ )x(c £ )Vc £ ■ Vn e | 

Jn 

+7 / |nJF £ '(n £ )x(c £ )Vc £ | + / nj|u £ |') ||v?||iyi.«>(n) 
for all t > 0 . Due to the embedding W 2,q (Q) f+ 1,00 (fl) for q > 3 , we deduce C\ > 0 such that 

J IK n e)lll(VE 2 >9(n))* ^ 

< Ci (^7(7-1)/ / D £ (n £ )nr 2 |Vn £ | 2 + 7 [ [ 71 £ (n, e )n 2 '~ 1 |Vra £ | 


0 Js] 


0 Jsi 


nj|u £ 


+7(7-1) / / |nj 1 i ? £ (n £ )x(c £ )Vc £ • Vn £ 

io in 

+7 [ [ \nfF' £ { n £ )x(c £ )Vc £ \ + / 

Jo Jn Jo Jn 

=: h + I2 + I3 + h + h 

for all T > 0. By (1.7), (5.9), (5.1) and Young’s inequality, we can find C 2 > 0 such that 
h < ^27(7-1)/ / (n £ + £) m n)T 1 |Vn £ | 2 


(5.17) 


< 


< 


^2(7- 1) 
m + 1 J 0 
712(7 - 1 ) 


/ 0 lo 

rT 



2 (m + 1 ) \J 0 jq 
< C 2 (T + 1 ). 


V(n e + e ) m+1 • Vn/ 


n e + e ) m+1 | 2 + 



I VnJ|‘ 


0 1!) 


(5.18) 
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Employing (1.7), (5.3) and Young’s inequality we have C3 > 0 such that 

rT r D £ (n E ) ~ ' r 


h < 


< 


UL 


'0 Jn 
< C 3 (T + 1), 

where we have used when 1 < rn < 2 

rT 


|Vn £ | 2 + 1 f [ D E (n E )np- ] 

n n e 2 y 0 Jn 

7 [ T f D e (n e ) l 

2./o Jn 


-|Vn e | 2 + 


7^2 f T 


(n e + e) m+27 


-2 


'0 in 


(5.19) 


/ f(n £ + e )"*^" 2 = [ f (n £ + e) m < C 4 (T + 1) by (5.2) 
Jo Ju Jo Jo 


with C4 > 0 and in the case m > 2 

rT 


Jo Jn 

From (2.6), (2.10), (4.2) and (5.1), we estimate 

j-T 

h < 


(n e + e) m+2l ~ 2 = f f (n £ + e) 2 ^" 1 ) < C(T) by (5.8). 
Jo Jn 


7(7 - 1 )C 5 [ [ \n] 1 

Jo Jn 


< Csi'y-l) / / | VnJ • Vc, 

Jo Jn 

C 5 ( 7-1) x ^ 


Vn E • Vcr 


< 


'o Jn 


|Vn * |2 + U 


< C e (T + 1) 


(5.20) 


with C5 > 0 and > 0. Applying (2.6), (2.10), (4.2) and (5.2), we find C 7 > 0 and C 8 > 0 such 
that 


h < C 7 ( f T f n/ + [ T [ |Vc £ 

\J 0 Jn Jo Jn 
< C 8 (T + 1). 

Moreover, we use (5.2) and (5.5) to give Cg > 0 and C]0 > 0 fulfilling 


(5.21) 


h < eg (r f n y + r [ 

\J 0 Jn Jo Jn 


\u £ \ 3 


/0 Jn 

< CioCT + l). 


(5.22) 


Then, (5.14) follows by combining (5.17)-(5.22). Multiplying the second equation in (2.1) by 2 ^~ 

with ip G C°°(i4) and the third by 0 € (Co^fi)) 3 , we obtain (5.15) and (5.16) in a completed 
similar manner (see [29] for details). □ 


6. Global weak solutions for (1.2)-(1.4) 

We are now in the position to construct global weak solutions for (1.2)-(1.4). Before going into 
details, let us first give the definition of weak solution. 

Definition 6.1. We call (n,c,u) a global weak solution of (1.2)-(1.4) if 

neLl c (nx[ 0,oo)), cG L^ c ([0,oo); W 1 ’ 1 ^)), u G (^([0, 00 ); Wq 1 ’ 1 ^))) 3 
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such that n > 0 and c > 0 a.e. in !1 x ( 0 , oo), and that 
™/(c) G Ll oc ([0, oo); L 1 (fl)), 

D(n)Vn, nx(c)Vc, nu and cu belong to [Lj oc ([0, oo); L 1 (fl ))) 3 , 
u®u£ {L} oc {[ 0,oo);L 1 (O))) 3x3 


and that 


roo r roo r 

/ / n t (j )i - nu- V0i = - 

Jo isi Jo Jn 


/* roo /* 

/ D(n)Vn ■ V0i + / / nx(c)Vc • V0i, 

JQ jo Jo 


roo /• 

7o 

/♦oo 


Q 02 - 


C« • V 02 = — 


/o Jf2 

roo 


roo r 

7o 


Vc • V 02 — 


[ nf{c)(j) 2 , 
Jn 


roo r roo r roo r roo r 

/ Uf(j) 3-K / / U®U-\74>3 = - / / Vu-V0 3 + / 

Jo id jo Jn jo id Jo Jn 


nV<f> • 03 


hold for all 0i € x [0, oo)), 02 G Co°(fl x [0,oo)) and 03 € (Co°(fi X [0,oo ))) 3 satisfying 

V ■ 03 = 0. 


We can now pass to the proof of our main result. 

Proof of Theorem 1.1. Let 7 be given by Lemma 5.1 and set 


0 : = 


^ 2 , 1 < m < 2, 


2 , 


m > 2. 


By Lemma 3.3, Lemma 5.1 and Lemma 5.3, for some C > 0 which is independent of e, we have 


— C r (^+ 1 )> 

( n e)tllLj oc ([0 ) oo);(Wa.9 ( n )) * ) ^ C ( T + X ) with SOme 9 > 3, 
V^llLf oc ([0,oo);tY 2 . 2 (n)) ^ C ( T + !)> 

(■v/^)t|| 5 5 < C(T + 1), 

' U £llL 2 oc ([ 0 ,°o);VE 1 . 2 (n)) - + 1 )> an ^ 


U Et | 


i 5 


iilc([ 0 ,oo);(w 0 ; ( ?(n))*) 


< C(T + 1) 


for all T > 0. Therefore, the Aubin-Lions lemma ( [12], see [16] for the case involving the space L p 
with p = 1 ) asserts that 

( n e)ee(o,i) is strongly precompact in L^ oc (£l x [ 0 ,oo)), 

(\/ce)ee(o,i) i s strongly precompact in Lf oc ([0,oo);W 1,2 (Q)) and 
( u e)ee(o, l) i s strongly precompact in L 2 oc (Cl x [ 0 ,oo)). 

This yields a subsequence e := £j G (0,1) (j G N) and the limit functions n, c and such that 

nj —> n 7 in L^ oc {pL x [ 0 ,oo)), and a.e. in x ( 0 ,oo), 

3m+2 _ 

n e ->• n in (fi x [ 0 ,oo)), 

3m+2 

D £ (n £ )Vn £ U(n)Vn in L™(ft x [0, oo)) 
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and 


\fcl \fc 


V C 4 


in L 2 oc ([ 0, oo); kk 1,2 (fl)) and a.e. in Q x (0,oo) 
in x (0, oo)), 

in Lf oc (ti x [0, oo)) 


as well as 


as e 


as e 


u £ —} u 
u £ ^u 

u £ 1 u 
Vu £ —*■ Vu 


in Lf oc (£l x [0, oo)) and a.e. in 17 x (0,oo), 

in L°°([0,oo);L^(n)), 

10 

in L i 3 oc ( n x [0,oo)), 
in Lf oc (n x [0, oo)) 

0. Moreover, using interpolation inequality for L p -norm, we have 

io 

n £ -A n in L^ c (17 x [0, oo)) 

0. According to (2.10) and the Lebesgue dominated convergence theorem, we obtain 
F e( n e)x{c £ )ci -A x(c)c* in L 2 o 0 c (fj X [0, oo)), 


/(c, 


■e; ^ /( c ) 
c £ c 


as e -A 0. Therefore, 


F £ (n £ ) - 

Y £ u £ - 

n £ u £ 


n 

u 


nu 


n £ F £ (n £ )x(c £ )Vc £ nx(c)Vc 
F e {n e )f{c e ) -A n/(c) 
c e u e -A eu 
Y £ u £ <8 u -A u <8 u 


hl ^oc( fi x [0, oo)), 

in Ll c (n x [0, oo)), 

10 

in L ioc^ x [0, oo)), 
in L 2 loc {Sl x [0, oo)) 

in x [0, oo)), 

in x [0,oo)), 
in x [0,oo)), 
in Lj oc (n x [0, oo)), 
in L} oc (p, x [0, oo)) 


as e -A 0. Based on the above convergence properties, we can pass to the limit in each term of 
weak formulation for (2.1) to construct a global weak solution of (1.2)-(1.4) and thereby completes 
the proof. 
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